Neutron scattering study of ferroelectric Sn2P2S6 under pressure 
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Ferroelectric phase transition in the semiconductor Sn2P2S6 single crystal has been studied by 
means of neutron scattering in the pressure-temperature range adjacent to the anticipated tricritical 
Lifshitz point {p ~ 0.18 GPa, T ~ 296 K). The observations reveal a direct ferroelectric-paraelectric 
phase transition in the whole investigated pressure range (0.18 — 0.6 GPa). These results are in a 
clear disagreement with phase diagrams assumed in numerous earlier works, according to which a 
hypothetical intermediate incommensurate phase extends over several or even tens of degrees in the 
0.5 GPa pressure range. Temperature dependence of the anisotropic quasielastic diffuse scattering 
suggests that polarization fluctuations present above Tc are strongly reduced in the ordered phase. 
Still, the temperature dependence of the (200) Bragg reflection intensity at p = 0.18 GPa can be 
remarkably well modeled assuming the order-parameter amplitude growth according to the power 
law with logarithmic corrections predicted for a uniaxial ferroelectric transition at the tricritical 
Lifshitz point. 
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I. INTRODUCTION 



Ferroelectric substances with a narrow electronic 
band gap have recently attracted considerable atten- 
tion due to their potentially interesting thermoelectricji 
photovoltaic,^''^ and other photoactive properties;^"— Uni- 
axial ferroelectric chalcogenides represent one of the 
best known families of ferroelectric semiconductors In 
particular, solid solutions of Sn2P2S6 and Sn2P2Se6 
have been recognized as extremely interesting model 
systemsi^"— All of them have an identical parent para- 
electric structure at high temperatures (P2i/n) and the 
same ferroelectric phase at low temperatures (Pn).— "— 
At ambient pressure, crystals with low Se/S concentra- 
tion ratio show a direct ferroelectric-paraelectric phase 
transition, while those of high Se/S concentration exhibit 
an intermediate incommensurate (IC) phase (see Fig.[T]). 
The IC phase evidenced in Se-rich solution o^'^" resem- 
bles a periodic array of anti-parallel ferroelectric layers 
with thickness of about 4nm {\qi\ « O.OSA"^). The 
lower-temperature phase transition in Sn2P2Se6 could be 
thus narrated as a discontinuous transformation from a 
strictly regular "ferroelectric nanodomain" arrangement 
to a structure with usual, macroscopic ferroelectric do- 
mains. This situation is known from many other so-called 
type-II IC ferroelectric systems, such as thiourea, BCCD, 
NaN02, etcJ^ The peculiarity of the Sn2P2(Se:^Si_2;)6 
solid solution is that here with decreasing value of Se/S 
concentration x, the temperature range of the interme- 
diate IC phase continuously decreases and finally the 
paraelectric-IC and IC-ferroelectric phase transition lines 
join near the so-called Lifshitz point (LP, see Fig.[T|) at 
X « 0.28 and T « 280 K. Theoretically, the modulation 
period should diverge therei^"— 



Later, it was predicted that a similar LP could also 
occur in the temperature-pressure (T — p) phase diagram 
of the pure Sn2P2S6 crystal. Several experiments indeed 
confirmed that the phase transition line splits into two 
lines after crossing the point with coordinates around 
p = 0.18 GPa, T = 296 K. The pure Sn2P2S6 system 
is obviously even more attractive than the x = 0.28 com- 
pound where, in principle, the LP-related phenomena 
could be obscured by compositional fluctuations inher- 
ent to solid solutions. In addition, thermodynamic anal- 
ysis of various physical quantities indicated that the Lan- 
dau free-energy expansion coefficient of the fourth-order 
term in polarization strongly decreases with pressure in 
Sn2P2S6f2 In fact, the thermodynamic coordinates of the 
tricritical point (TCP), at which this quadratic term van- 
ishes, were predicted to fall in the vicinity of the above 
discussed LP. This suggests that the fortunate intersec- 
tion of two special lines - the line of the TCPs and the 
line of LPs, known also as a tricritical Lifshitz point 
(TCLP) - might be reached in Sn2P2S6 under hydrostatic 
pressure.— 

However, several conflicting statements were reported 
about the nature of the high-pressure phase transition 
of Sn2P2S6. In particular, recent analysis of the speciflc 
heat anomalies^- supports the earlier guesa^i that the 
conditions for the LP are met at lower pressure than the 
conditions for the TCP. This means that only the LP 
is realized in the pressure-temperature phase diagram 
of Sn2P2S6. This scenario was repeatedly assumed in 
Rcfs..28,-32 and implies that the special point with coor- 
dinates around p ~ 0.18 GPa, T ~ 296 K is of the same 
nature as the LP in the temperature - concentration dia- 
gram of Sn2P2(Sea;Si_a;)6 solid solution (see Fig.[T]). On 
the contrary, other recent measurements, such as bire- 
fringence measurements of Ref.H^ were used to support 
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FIG. 1: (Color online) Schematic temperature-composition 
phase diagram for the system Sn2P2(Sei:Si-i)6 according to 
Refs.[9j,18,. Dashed line - second-order phase transition; solid 
line - first-order phase transition. 



(a) p = 0.6 GPa, 7= 220 K (b) p = GPa, 7= 383 K 
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FIG. 2: (Color online) Accessible part of the a*c* scatter- 
ing plane with schematic representation of neutron scattering 
intensities of principal Bragg reflections in the paraelectric 
phase. The areas of the circles are proportional to abso- 
lute values of neutron scattering structure factors (a) in the 
present experiment at p = 0.6 GPa and T — 220 K and (b) 
p = GPa and T — 383 K values recalculated from the x-ray 
data of Ref.fsJ using the Jana2006 computer program^ with 
the neutron wave number set to kf = 1.48 A~^. 



the original scenario of the canonical theory of Ref.ls^: 
That the observed special point in the T ~ p diagram of 
Sn2P2S6 is merely a TCP but surely not a LP.^^ More- 
over, the absence of the LP in the T — p phase diagram 
of Sn2P2S6 was also inferred from the comparisons with 
the Blume-Emery-Grifhth's modeljS^ On the top of that, 
there is a considerable disagreement about the domain 
of stability of the high-pressure intermediate phase of 
Sn2P2S6 among various sources.— i^i^i^ 

Despite these challenging issues, diffraction investiga- 
tions of the intermediate phase of Sn2P2S6 were not re- 
ported so far and there is no direct evidence of its IC 
nature. The aim of this work is to bring at least a partial 
answer to the above raised issues by discussion of new 
results from high-pressure investigations of Sn2P2S6 by 
means of neutron scattering. 

II. EXPERIMENTAL DETAILS 

The neutron scattering experiment was performed at 
the Laboratoire Leon Brillouin (Saclay, France) on the 
2T1 thermal neutron beam and 4F2 cold neutron beam 
three-axis spectrometers (TASs). Spectrometers were op- 
erated with PG(002) analyzer set to the neutron wave 
number kj = 2.662A~^ and kf = 1.48A~-'^, respectively. 
The energy resolution at zero-energy transfer was 0.2 and 
0.05 THz, the Q resolution about 0.01 A'^ and 0.005 A'^ 
on the 2T1 and 4F2 TASs, respectively. 

The Sn2P2S6 single crystal used in this experiment was 
grown from melt by the Bridgeman method at Uzhgorod 
University. The sample (w 0.8 cm'^) was placed in a spe- 
cial helium pressure cell designed for neutron scattering 
experiments in a temperature range of 20 — 305 K and 
under hydrostatic pressures up to about 0.6 GPa. 

The crystal was mounted so that the scattering plane 
matched its crystallographic mirror plane niy (through- 



out the paper, we stick to the pseudo-orthorhombic unit 
cell defined by a = 9.378 A^= 7.488 A, c = 6.513 A, 
and (3 = 91.15°, as in Ref.lSa)- This horizontal plane 
thus also contains ferroelectric polarization direction (at 
about 10 — 15° out of the a axis). A reasonably ac- 
cessible part of the reciprocal space in such a scatter- 
ing geometry is shown in Fig. [21 A rough survey of un- 
corrected neutron diffraction intensities of main Bragg 
reflections in Sn2P2S6, as detected at T = 220 K and 
p — 0.6 GPa in our experiment, is indicated by variably 
sized symbols in Fig.t^fa). Results agree fairly well with 
the neutron diffraction intensities calculated from struc- 
tural data of the ambient-pressure paraelectric phase of 
Sn2P2S6 [Fig.Hb), data at T = 383K of Ref.ig. 

The a*c* plane of Fig. [2] is the same a*c* plane within 
which the satellite reflections with a modulation wave 
vector Qi « 0.085c* — O.Ola*, demonstrating the IC 
phase, were observed in the Sn2P2Se6 crystal.^'' Indepen- 
dent x-ray and neutron scattering experiment c'^^'^^ agree 
that satellite reflections of Sn2P2Se6 are particularly 
strong for large momentum transfers Q, roughly parallel 
to the spontaneous polarization (Q || Pg). This is a con- 
sequence of the usual situation caused by the transverse 
character of the frozen polarization wave (qi_LPs), and 
the same can be anticipated for the supposed Sn2P2S6 
satellite reflections. Clearly, detection of satellites with a 
very small modulation wave vector, comparable to the 
tails of the Q resolution of the instrument, could be 
problematic in the vicinity of a strong Bragg reflection. 
Unfortunately, systematic absences are common to both 
main and satellite reflections, that is, the satellites are 
expected only around the allowed reflections in the a*c* 
plane. For this reason, most of our measurements were 
taken in the vicinity of the (200) Bragg reflection. This 
is an allowed reflection, but rather a weak one, and at 
the same time the corresponding scattering vector Q is 
almost parallel to Pg. 



3 



a. 
>- 

'U3 
C 



■ I—' 



1 

-I—' ' 



(a) 



p = 0.18 GPa 
heating, 2T1 




□ data 

power law (p 



250 260 270 280 290 
Temperature (K) 



300 310 





4- 


CD 
Q. 

>. 




tensit 


3- 








2- 


>< 




tensi 


1 - 


c 





(b) 



p = 0.6 GPa 
heating, 2T1 




□ data 

power law (P = 0.160) 





160 170 



180 190 200 
Temperature (K) 



210 220 



FIG. 3: (Color online) Temperature dependence of the (200) 
reflection intensity (a) at O.fSGPa and (b) at 0.6 GPa. The 
full points and solid lines refer to the fits using Eqs. @ and 
([5]) with and without logarithmic corrections, respectively. 



III. RESULTS 

In agreement with our anticipation, the neutron 
diffraction intensity of the (200) Bragg reflection in the 
paraelectric phase was several times smaller than in the 
ferroelectric phase at ambient conditions. In fact, mon- 
itoring of the (200) reflection allowed us to probe quite 
accurately the ferroelectric order parameter. As an ex- 
ample, the temperature dependence of the (200) Bragg 
reflection intensity at 0.18 GPa and 0.6 GPa is shown in 
Figs.|3^a) andin^b). The intensity visibly starts to rise at 
a well defined temperature that can be associated with 
the ferroelectric phase transition. The temperatures thus 
obtained are given in Table I and within the possible ex- 
perimental error in the pressure and temperature deter- 
mination they correspond well to those derived from the 
compressibility anomalies reported in Ref.[28l. 

In the vicinity of these transition temperatures, we 
have searched for possible indications of a satellite reflec- 
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FIG. 4: (Color online) Zero-energy-transfer scans in the (200) 
Brillouin zone of the Sn2P2S6 crystal taken in the vicin- 
ity of the ferroelectric phase transition (at 0.6 GPa; near 
Tc = 205.3 K). (a) Scans in the [-0.03,0,1] direction (ap- 
proximately along the diffuse streak), and (b) scans in the 
[1, 0, 0.05] direction (across the diffuse streak). Measured scan 
trajectories are schematically indicated in the insets. Solid 
lines are just guides for the eye. 



tion. The pressure of 0.18 GPa is supposedly that of the 
TCLP, while the measurements at 0.3 GPa and 0.6 GPa 
are taken in the region where finite temperature range 
of the IC phase was anticipated. For example, we have 
recorded intensity maps in a 2D area around the (200) 
Bragg refiection, extending up to ±0.2 r.l.u. along both 
a* and c* at 0.6 GPa and temperatures of 205.5, 206, and 
206.5 K (at 0.2, 0.7, and 1.2 K above Tc = 205.3 K) and 
similar maps were recorded also at 0.3 GPa and temper- 
atures 267, 270, 273, and 276 K (Tc = 271 K). However, 
no obvious satellite peaks have been seen. 
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FIG. 5: (Color online) Schema of reciprocal-space trajectories 
of data scans shown in Fig.|4l superposed with contour plots 
of the critical diffuse scattering known from ambient pressure 
x-ray diffraction measurements of Ref.l4ll. Full triangles in- 
dicate positions of incommensurate satellite reflections of the 
ambient-pressure incommensurate phase of the related com- 
pound Sn2P2Se6, as reported in Ref.[20l. 



At the same time, our observations disclosed traces of 
weak and broad diffuse scattering ridges, extending over 
the "soft" direction, roughly perpendicular to the direc- 
tion along which the spontaneous polarization is formed. 
FigureSlJa) shows an example of a scan taken along this 
diffuse scattering ridge at several temperatures. Exam- 
ples of transverse scans are shown in Fig.H^b). Trajecto- 
ries of these scans are schematically shown in Fig. [3 along 
with the contour plots of the ambie nt- p ressure diffuse 
scattering topography taken from Ref.|4ll Note that the 
diffuse scattering ridge corresponds to the smaller peak 
in Fig.mjb), while the larger intensity peak is merely a 
spurious leakage of (200) Bragg reflection scattering due 
to the strongly anisotropic tails of the instrumental res- 
olution function. 

Some of the scans displayed in Fig.|H^a) show a non- 
monotonic intensity decay with the distance from the 
Brillouin zone center. However, due to the consider- 
able background intensity caused by high-pressure-cell 
environment, this can hardly be considered as a signif- 
icant indication of incommensurate satellite reflections. 
Rather, the anisotropy of this diffuse scattering is ob- 
vious and strikingly similar to the critical fluctuations 
previously seen by x-ray scattering at ambient pressure 
conditions 4i The assignment of the diffuse scattering to 
the order-parameter fluctuations is further corroborated 
by the observed increase in its intensity as the transi- 
tion is approached from above and by its rather abrupt 
reduction just below the phase transition point [see full 
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FIG. 6: (Color online) Temperature dependence of both elas- 
tic and quasielastic scattering in the (200) Brillouin zone at 
0.6 GPa: 

(a) selected constant-Q scans for several temperatures, (b) 
principal Bragg reflection intensity at q = (0, 0, 0) and elastic 
diffuse scattering at q = (0.02, 0, -0.07). 



points in Fig.[SJb)]. Moreover, few constant-Q spectra 
taken at selected positions in the high-temperature phase 
[Fig-inia)] suggest that this diffuse scattering is quasielas- 
tic in its nature, with FWHM of the order of 0.2 THz. It 
will be shown in the next section that several impor- 
tant conclusions can be drawn from such measurements, 
even though the high elastic and inelastic scattering back- 
ground did not allow their more quantitative analysis. 



IV. DISCUSSION 

Absence of the TCLP, three alternative scenarios. 
Summarizing, our results do not provide any support for 
the LP in the T — p phase diagram of Sn2P2S6 up to the 
hydrostatic pressure values of about 0.6 GPa. Thus, the 
most natural conclusion seems to be that there is no IC 
phase in this T — p region at all (" scenario A" ) . The 
absence of any clear evidence for the LP in the above 
displayed observations led us nevertheless to perform de- 
tailed and careful comparisons of our results with those 
of previously reported experiments, which did confirm 
the presence of the LP. In particular, we have traced 
the T — p phase diagram of Sn2P2S6 (see Fig.[71) accord- 
ing to the reported anomalies in the compressibility data 
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FIG. 7: (Color online) Temperature-pressure phase diagram 
of Sn2P2(Se^Si-^)6 solid solution.^^i^i^i^i^ The shaded ar- 
eas indicate intermediate phases between high-temperature 
paraelectric and low-temperature ferroelectric phases for 
given concentrations x. The stars stand for Lifshitz points 
of solid solution with x equals: 1, 0.28, 0.2 0.04, and (from 
left to right). Temperature ranges of neutron scattering mea- 
surements reported in this work are indicated by vertical rect- 
angular areas. 



(Ref.lH, circles in Fig.[71) and in the dielectric permittiv- 
ity data (Ref.H^, triangles in Fig. [7]). Both these measure- 
ments suggest splitting of the phase transition in Sn2P2S6 
{x = 0) at around p = 0.18 GPa, T = 296 K, which was 
naturally interpreted as the expected LP. At the same 
time, the stability range of the anticipated intermedi- 
ate phase is markedly different in these two measure- 
ments: It opens at the rate of about 190 K/ GPa accord- 
ing to Ref.(43l ("scenario B"), while there is only about 
20 K /GPa opening rate deduced from measurements of 
Ref.lii ("scenario C"). 

It is also instructive to compare these data with the 
T — p phase diagram of the Sn2P2Se6 crystal, in which 
the IC phase was firmly established at least at ambi- 
ent pressure. In fact, the temperature range of the IC 
phase of Sn2P2Se6 (limited by pairs of square symbols in 
Fig.[71) opens at a rate of about lOOK/GPa, which is just 
an intermediate value between the two visibly conflicting 
values given for Sn2P2S6 (20 and 190 K/GPa, respec- 
tively). Let us note that the similar data for the x = 0.2 
solid solution crystals (cross symbols. Ref. i27h and the 
X = 0.05 solid solution crystals (Ref.l27l not included in 
Fig. [7]), suggest already a quite narrow range of the IC 
phase, which seems to be much more consistent with the 
narrow IC phase for Sn2P2S6 (scenario C). The latter 
tendency could be also easily associated with the fact 
that the IC phase temperature range width is strongly 
decreasing with the S/Se ratio at ambient pressure (see 
Fig.[T]). Nevertheless, the phase stability of the IC phase 



is a strongly nonlinear function of the Se concentration 
at ambient pressure (see Fig.[T]) and so it could be a non- 
monotonous function of the Se concentration at higher 
pressures. 

In spite of these disagreements about the tem per ature 
range of the IC phase in Sn2P2S6, results of Ref.[28l (sce- 
nario C) and Refs. fislH (scenario B) agree well about 
its onset at about 0.18 GPa, T = 296 K. This location of 
the LP appears to be fairly well consistent with T — p 
coordinates of LPs reported for the several solid solu- 
tions with different Se concentrations (LPs are denoted 
by star symbols in Fig. [7]). Likewise, the "theoretical" LP 
of the Sn2P2Se6 compound, obtained by extrapolating of 
the experimental paraelectric-IC and the IC-ferroelectic 
phase transition temperatures towards negative pressure 
values (see Fig. [7]), seems to lie on a more or less straight 
line in the projected T — p diagram. 

Elimination of the scenario B by continuity arguments. 
So how credible is the outcome of the present result with 
respect to the previous experiments? Perhaps the most 
significant advantage of the present experiment is that 
here we employed a scattering technique, which allows 
us to appreciate the amplitude of the order parameter as 
well as its wave vector. In particular, it is important that 
our experiments positively prove that the zone-center 
ferroelectric order parameter vanishes at temperatures 
given in Table I in a fairly continuous manner. It is well 
known that virtually all displacement-type IC dielectrics 
obey the standard continuity paradigm supported by the 
Landau theoryiS: (i) the phase transition from the para- 
electric to the IC phase is a continuous transition, and 
(ii) the modulated- ferroelectric "lock-in" transition is a 
first-order phase transition but with a relatively mod- 
erate jump of the modulation amplitude. Thus, already 
from these results, there is no space for an IC phase above 
temperatures given in Table I. Based on this argument, 
we can safely refute the phase diagram with a broadly 
opened IC phase^i^ (scenario B). 

In order to make our conclusion even more convincing, 
let us recall that the atomic positions Vj in the vicinity 
of the ferroelectric structural transition can be approxi- 
mately expressed through the paraelectric reference ones 
(rj,para) by means of the frozen order-parameter ampli- 
tude 77 and corresponding eigenvector ej as 

Tj — Tj.para -I" V^j ■ (1) 

Inserting this into the standard expression for the elastic 
neutron scattering structure factor 

F{Q)=Y,bje'^-''^e-^' , (2) 

j 

where bj is the neutron scattering length, Wj is the 
Debye- Waller factor of the jth atom, and Q is the scatter- 
ing vector and keeping the leading term in rj, one obtains 
the reliable and frequently exploited relation 

F(r;,Q) «F(0,Q)(l + iryQ-ej). (3) 
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This expression implies that the order parameter is mod- 
ifying the elastic coherent neutron scattering intensity as 

= /(0,Q) + 772/(0, Q)|Q.e,f . (4) 

Thus, the extraneous contribution to the Bragg reflection 
intensity is proportional to the square of the ferroelec- 
tric order-parameter amplitude (77^). Moreover, it can 
be shown, that for a long-wavelength IC phase of type 
II, the second term in the above expression also quite 
well approximates the satellite reflection intensity. Con- 
sequently, the intensity by which the principal reflection 
is increased just below the lock-in transition should be 
comparable to the intensity of the satellite reflection just 
above the lock-in transition. 

In other words, one may expect that above the lock-in 
phase transition, the satellite reflection in the vicinity of 
the (200) Bragg reflection should have about the same in- 
tensity as the order-parameter-related part of the (200) 
Bragg reflection intensity itself. Here in our measure- 
ments (see Fig. [3]) the anomalous intensity of the (200) 
Bragg reflection decreases considerably and continuously 
as the temperature approaches the transition tempera- 
ture (given in Table I), from below, and thus we know 
that the order parameter vanishes there completely. This 
suggests that such transition cannot be a lock-in transi- 
tion, and a modulated phase persisting above this tem- 
perature as proposed from the diagram of Refs.H^H^ 
(scenario B) can hardly be expected. 

Comparison of scenarios of A and C. Therefore, the 
only possibility on how to reconcile our observation with 
a hypothesis of an IC phase is to assume that it ex- 
ists below the transition temperatures of Table I. This 
could happen for a modulation wave vector smaller than 
the momentum resolution, as then the satellite reflec- 
tion intensity would be integrated into the (200) Bragg 
reflection scattering. Since a systematic broadening of 
the (200) Bragg reflection was not noticed in our exper- 
iments, it would mean that the modulation period is of 
the order of 100 nm. The phase transition temperature 
of Table I would correspond to a paraelectric-IC phase 
transition Tj. This is in an acceptable agreement with 
phase transition temperatures given by Ref.[2^ (scenario 
C). 

From the same woik^ one may infer that the lock- 
in transitions Tiock-in would fall at about 3 K (5 K) be- 
low Ti at 0.3 GPa (at 0.6 GPa). Obviously, a small jump 
of the order-parameter amplitude at these temperatures 
can hardly be identifled in our data. In principle, this 
scenario seems to be compatible also with the birefrin- 
gence data in Ref. 45. There are indeed tiny jumps in 
the temperature dependence of the anomalous part of the 
birefringence of Sn2P2S6 at high pressures, even though 
the data alone were not considered as significant.^^ At 
the same time, very similar jumps were observed at the 
well established Tiock-in temperature in the temperature 
dependence of the birefringence data of Sn2P2Se6, and 
they are also quite smallj^"— despite a much broader 



TABLE I: Phase transition temperatures and /? parameters 
of formulas (U) O without any logarithmic correction for the 
(200) Bragg reflection measured at given pressure p in the 
ferroelectric phase on the 2T1 TAS. 

p (GPa) Tc (K) 13 heating rate (K/min) 

0.18 299.6 0.215 0.22 

0.3 271.0 0.196 0.06 

0.6 205.3 0.160 0.22 



TABLE IL List of critical dimensions du and critical expo- 
nents for uniaxial ferroelectrics according to Refs. [22 23. The 
star symbol * indicates the logarithmic corrections in the 
power laws. (U-uniaxial dipolar, L - Lifshitz, T - tricritical, 
m-dimensions of Lifshitz subspace). 

system du ce j3 y x 

U 3 0* I* 1* i 

rp o liH^l^tii^f: 1 

'- "'24^3 
UT 2i i i 1 - 

ULT, m = 1 3 ^* J* I* Jo 



range of the IC phase there and presumably more abrupt 
lock-in transition. Let us note that in contrast, the com- 
pressibility data of Ref.[2^ show a larger anomaly at the 
presumed Tlock-in transition, while Tj is barely notice- 
able. The compressibility anomalies, however, do not 
directly scale with the order-parameter amplitude. On 
the contrary, it is likely that the compressibility anoma- 
lies are sensitive to the modulation period, similarly as it 
was demonstrated for example in the case of modulated 
phases of BCCD.'*^ Therefore, the large compressibility 
anomaly at Tjock-in might be compatible with scenario 
C as well. 

Several possible experimental problems that might be 
at the origin of the contradictory reports on the splitting 
of the phase transition in Sn2P2S6 have been recently 
discussed for example in Ref. 45- Perhaps the most seri- 
ous is the peculiar sensitivity of Sn2P2(Sea;Si_a;)6 solid 
solutions to the absorption of light. It was reported 
that an a; = 0.2 crystal irradiated by 0.5 W krypton 
(647 nm) laser radiation shows two close anomalies in the 
temperature dependence of the light transmission even 
though the crystal with this composition has nominally 
a single continuous phase transition. We speculate that 
this partly understood photovoltaic or photoelastic phe- 
nomenon could be responsible for the extraneous anoma- 
lies that are behind "scenario B" (experiments reported 
in Ref.i43i). However, strictly speaking, same worries ap- 
ply to the experiments supporting the narrow IC phase 
(scenario C). 

Order- Parameter fluctuations. Additional support for 
scenario A can be derived from the observations of the 
temperature dependence of the critical fluctuations of 
the order parameter. In displacive ferroelectrics, the in- 
tensity of the associated quasielastic diffuse scattering is 
typically proportional to l/wgjyj, where wsm is the soft 
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mode frequency at a given wavevector. Therefore, the 
increase of the critical scattering in the vicinity of the 
phase transition is directly associated with soft branch 
frequency lowering. In the incommensurate phase, the 
soft mode branch is split into two branches, called pha- 
son branch and amplitudon branch, respectively.—!^ Ob- 
viously, the diffuse scattering is mainly determined by the 
lower frequency phason branch. ^'^ Since for small modu- 
lation amplitudes the phason gap remains zero, the pha- 
son branch and the diffuse scattering does not show much 
temperature evolution within the range of incommensu- 
rate phase liSiSiH 

The critical quasielastic diffuse scattering observed 
within a few degrees above the phase transition in 
present experiments (at 0.6 GPa, see Figs.H] andlH]) can 
be indeed well understood as due to an overdamped 
soft mode branch. In particular, the pronounced tem- 
perature dependence of the diffuse scattering intensity 
shown in Fig.|6l^b) (data corresponding to the right-hand 
side scale) suggests soft mode frequency lowering upon 
approaching the nominal transition temperature from 
above. Such temperature dependence is typical for a 
paraelectric phase, but not for an incommensurate phase. 
At the same time, the steep intensity drop below Tq sug- 
gests marked stiffening of the soft mode frequency. This 
can be indeed expected for a ferroelectric phase transition 
with a partial order-disorder character (see discussion in 
Refs.mii) and it is in agreement with direct measure- 
ments of both the static dielectric permittivity temper- 
ature dependence and the soft-mode frequency temper- 
ature dependence (Ref.lSTf) near ambient-pressure phase 
transition of Sn2P2S6. 

As already mentioned, the quasielastic, that is, dynam- 
ical nature of this diffuse scattering intensity is apparent 
from Fig-inja). Spectra taken in the immediate vicinity 
of Tc suggest a bell-shaped quasi-elastic intensity compo- 
nent with FWHM of the order of 0.2 THz, that is, visibly 
exceeding the instrumental energy resolution. Dynamic 
critical fluctuations with relaxational frequencies of the 
order of 0.2 THz would fit quite well with typical relax- 
ation frequencies expected for soft mode fluctuations in 
the system (for an overdamped soft phonon mode, the 
relaxational frequency is about ^^sm/ism, where 7sm is 
the damping of the soft phonon mode.) Moreover, the 
quasielastic signal which for temperatures near Tq would 
reach the values of the order of 1000 counts at zero en- 
ergy transfer in Fig.|6lja) (extrapolation of full symbols 
to zero-energy) does match quite nicely the 200 counts 
per minute of the q-dependent and T-dependent diffuse 
signal apparent from Figs.S] andlHl^b). Thus, in brief, all 
the observations considered here are in a rather good 
agreement with the scenario of a simple paraelectric- fer- 
roelectric phase transition (scenario A). 

Critical behavior of the order parameter. There has 
been considerable interest in theoretical understanding 
of the critical behavior of uniaxial ferroelectrics near the 
TCLP and many previous studies analyzed phase tran- 
sition anomalies of Sn2P2S6 and derived solid solutions 



with the hope to determine the critical exponents rele- 
vant to its unconventional universality class. Since the 
T — p phase diagram itself is far from being established, 
we did not try to enter into the problem of the criti- 
cal phenomena. Nevertheless, for the sake of comparison 
with other studies, we have analyzed our order-parameter 
data with the usually assumed power law dependence on 
the reduced temperature r = (Tc ~ T)/Tc. Interest- 
ingly, the critical exponent (3 obtained from the adjust- 
ment of the formula ^ with 77 oc r'^ to the experimental 
data (see Fig. [3]) was substantially lower than the mean- 
field value of 0.5. This result is in contradiction with 
the analysis of Refs.[2^ and [s^- On the other hand, our 
values of /3 measured above ptcp = 0.18 GPa are nicely 
joining those obtained below ptcp from the spontaneous 
polarization (respectively pyroelectric coefficient) data of 
Sn2P2(Seo.o4So.96)6-— Moreover, by far the best fit [see 
Fig.[3fa)] was obtained for adjustment to the expected 
logarithm-corrected power law^i^i^ 

r7 = AT^|ln(r)r, (5) 

with /3 — 1/4 and x — 1/10, that is, values imposed 
according to the predictions for the critical behavior in 
the vicinity of a TCLP of an uniaxial ferroelectrio^Si (see 
Table II.). 

V. CONCLUSION 

In conclusion, this paper clarifies the nature of the fer- 
roelectric phase transition in the Sn2P2S6 single crystal 
under a moderate hydrostatic pressure (up to 0.6 GPa). 
The neutron diffraction experiments indicate clearly that 
the ferroelectric order parameter vanishes in a fairly con- 
tinuous manner at a well-resolved transition tempera- 
ture. In spite of considerable efforts, no evidence of an in- 
termediate incommensurate phase was obtained. While 
we cannot fully exclude the existence of an extremely 
long-wavelength modulation below the phase transition 
temperature, the most simple explanation of our obser- 
vations is that at least up to 0.6 GPa, there is only a sim- 
ple paraelectric-ferroelectric phase transition in Sn2P2S6. 
Still, the temperature dependence of the order param- 
eter happens to be best reproduced by the logarithm- 
corrected power law expected for a uniaxial ferroelectric 
system in the vicinity of a TCLP. Thus, it is quite possi- 
ble that the thermodynamical properties of the Sn2P2S6 
compound and derived solid solutions may reflect the the- 
oretical closeness to the ideal TCLP behavior and so it 
should be taken into account in analysis of the data. 
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